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The nature of the ground states for a system composed of two coupled cavities with each con- 
taining a pair of dipole-coupled two-level atoms are studied over a wide range of detunings and 
dipole coupling strengths. The cases for three limits of exact resonance, large positive and negative 
detunings are discussed, and four types of the ground states are revealed. Then the phase diagrams 
of the ground state are plotted by choosing three different "order parameters". We find that the 
phase space, determined by the combinative action of detuning and the dipole coupling strength, is 
divided into four regions. This is different from the general Bose-Hubbard model and more richer 
physics are presented in the two-site coupled cavities system. That is, the insulator region may be 
polaritonic or atomic and the superfluid region may be polaritonic or photonic in nature. 

PACS numbers: 42.60.Da, 03.65.Yz, 42.70.Qs, 71.15.Ap 



INTRODUCTION 



The simulation of the strongly correlated many-body 
systems described by Bose-Hubbard model has received 
great advances in optical lattices [Tj-Hoj and coupled- 
cavity systems pj]-[2^. Both of them depend on the com- 
petition between the local interaction and the nonlocal 
tunneling, but there are also some differences for this 
two basic models. In optical lattices, the quantum phase 
transitions (QPT) in a gas of ultracold atoms with peri- 
odic potentials is described by the Bose-Hubbard model 
with on site two atoms interacting and hopping between 
the adjacent sites. However, in the coupled-cavity sys- 
tems, two types of particles are involved to study the 
many-body dynamics and its realization relays on the 
strong light-matter coupling regime. Therefore, the QPT 
is due to the transferring of the excitations from polari- 
tonic to photonic rather than purely bosonic or purely 
fermionic entities. However, the realization of the strong 
coupling in experiment is the greatest bottleneck for the 
QPT manipulated in coupled-cavity system. To be opti- 
mistic, with the progress in the realization of strong light- 
matter coupling regime in both atomic [25|, |26j and solid- 
state [27], [HI cavity quantum electrodynamics devices 
with single two- level emitters in high-Q resonators, the 
QPT for coupled- cavity arrays of Jaynes-Cummings(JC) 
model systems and its variants attracts more and more 
attentions. 

Subsequent works that deal with coupled nonlinear 
cavities arrays have addressed the dynamics in the two 
coupled cavities for its great freedom and flexibility, 
which can provides a convenience controllable platform 
for engineering the transport of quantum states via pho- 
tonic processes and the exact numerical solutions can be 
easily found and then some analytical approximations 
can be used. Only very recently the investigation of such 
systems of two coupled cavities has begun. The quan- 



tum states transfers [23|, the atomic state transfer [30], 
the quantum phase gates [3lj, the bipartite entanglement 
entropies j32| , the one-excitation dynamics [33|, the pho- 
ton correlatio ns 153 , the time evolution of the popula- 
tion imbalance |35j, the photonic tunneling effect [36j and 
the emission characteristics [37J have been studied. To 
the best of our knowledge, most of the previous works 
focused on this system are limited to the single atom- 
cavity interactions without consideration of an additional 
interatomic coupling. With the advances of the technol- 
ogy, the interatomic dipole-dipole interaction may be re- 
alized in several solid-state systems, such as an ensemble 
of quantum dots[38| and Bose-Einstein condensates [39| . 
At present, the stationary entanglement (iol Elj and 
QPT [42] for dipole-coupled two-level atoms in single- 
mode cavity have been investigated theoretically. Then 
an investigation of the QPT to the system of two coupled 
cavities with dipole-dipole interaction is also of consider- 
able significance and highly called for. 

The previous work[29| has identified some of the unique 
features of QPT in the coupled two-site JC model. The 
nature of the ground states for the system can be divided 
into four types corresponding to different parameter val- 
ues of atom-field detunings and cavity-cavity hopping 
strengths. Differing from the Bose-Hubbard model, the 
insulator state may be either atomic or polaritonic, while 
the superfluid state may be photonic or polaritonic in na- 
ture. In this paper, we extend the work[29j to each cav- 
ity containing two dipole-coupled atoms. We find that in 
the presence of weak cavity-cavity coupling, the effective 
on-site repulsion not only attributes to the atom-photon 
interaction and the detuning, but also depends on the in- 
teratomic dipole-coupled strength. So the atomic dipole- 
dipole interaction provides an additional parameter, and 
more importantly, richer physics for the system. 

The paper is organized as follows. We first describe 
the model under consideration and then simplify it to an 
effective form, with reformed atomic energy and atom- 
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field coupling strength. We then address the nature of 
the ground state of the system under a wide range of de- 
tuning and dipole coupling strengths. Thirdly, we discuss 
the quantum phase transitions of the system by choosing 
three different "order parameters". We finally present 
our conclusions. 



II. MODEL 

The system under consideration consists of two identi- 
cal single-mode cavities with each cavity containing two 
coupled two-level atoms through dipole-dipole interac- 
tion. The two cavities are coupled by hopping strength 
A, therefore the photons may hop between them. Be- 
sides, the model is ideal without taking into account the 
dissipation induced by atomic spontaneous emission and 
photonic escape from the cavities. In such a case, the 
Hamiltonian for the coupled two-cavity system is given 
by (h = 1) 



H 

Hi 

Ho 



H 
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Hi + H 2 + H\2-> 

uj c a\ai + i^aCrjcTi + g(a[ai + a\o\) 

i=l,2 

J(a\a 2 +oic4), 

3=3,4 
J(ajcr 4 + cr 3 cr|), 

A(a\a 2 + aia\), 



(2.1) 



where uj a and uo c are the resonance frequencies for atoms 
and cavities, respectively. Supposing the same atom-field 
coupling strength in the system, it is defined by an uni- 
form parameter g. a\ and are the creation and annihi- 
lation operators of the field in cavity i (i = 1, 2). cr+ and 
<jj represent the atomic raising and lowering operators of 
the atom j (j = 1, 2, 3, 4). For static atoms, the coherent 
dipole-dipole interaction between them can be given by 



J = |d| 2 (l-3cos 2 6)/? 
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(2.2) 



where, 1*12=1*1 — v 2 is the distance between the two atoms 
located at 1*1 and r 2 . is the angle between 1*12 and the 
atomic dipole moment d. Here, we assume the dipole 
moments of the two atoms are parallel to each other and 
are polarized in the direction perpendicular to the inter- 
atomic axis. Then, J can be simplified as 



J 



|d| 2 /r? 2 , 



and its strength can be adjusted by changing the posi- 
tions of the two atoms in each cavity. Hi and H 2 show 
the atom-field interaction and the atom-atom coupling in 
each site, respectively. The cavity-cavity coupling is de- 
picted by Hi 2 . The total excitation for the Hamiltonian 
H can be defined as TV 



assume the total number of excitations TV is conserved 
and exactly two excitations in the system. Then Hamil- 
tonian Hi and H 2 can be transformed into two simple 
forms by unitary transformation [40] , 

H leff = U{HiUi = u c a\ai + (u a + J)cr\(Ti 

+ V2g(a\ai + a\o{)], 

H 2e ff = U\HiU 2 = uj c a\a 2 + (u a + J)crja 3 

+ V2g(4a 3 + a 2 ai)}. (2.4) 



where, Ui exp[— ^(cj\a 2 -\-a\ai)\, U 2 exp[-|(crjcr 4 + 

04(73)]. In the transformed form, the dipole coupled 
atoms are denoted by two fictitious atoms. Only one of 
them couples to the field mode with frequencies uo a + J, 
but the other atom freely evolves decoupling from the 
field. The effective coupling strength also changes from 
g to \plg. In this paper, we pay our attention to the 
strong atom-cavity coupling regime which can be put into 
practice only when A <C g. On this condition, the eigen- 
states of the individual cavity should be expressed by the 
dressed states 



|0i 



\9i)\0i 



\ n i) = sin yl e i)l(™ - l)i> - cos y \9i)\rn), 



\nf) = cos y le^Kn- l)i> + sin y 1^)1^). (2.5) 

where i = 1, 2 indicates the cavity number, y/n is a pho- 
ton number state, n = a,rct8Lr^2^/2g^/n/(A + J)i ano ^ 
A = uj a — uj c is the detuning between the atom and the 
field. The eigenenergies of these eigenstates are 



E° = 0, 
= nt 



A + J 1 



T-V(A + J) 2 + 8 5 2„.(2.6) 



The effective form of Hamiltonian H can be rewritten 



H eff = H leff + H 2eff + H 12 



(2.7) 



— a\ai 



-a\a 2 - 



<j\<7i + g\g 2 . We 



Because there are only two excitations in the system, 
we can write out the the corresponding state of H e ff in 
the order of increasing energy and divide them into five 
groups, defined as |0i), |</> 2 ), Ifo), ^4), ^5), correspond- 
(2.3) ing to subspaces {|1^) <g> |1^>} ? {|2^) <g> |0 2 ), |0i) <g> |2^ )}, 
® |lj>,|lj"> ® |la », {|2j"> ® |0 2 >,|0i> |2+», 
{\li) II2")}, respectively. Obviously, the energy dif- 
ference between the adjacent subspaces depends on the 
parameter values in H e f f . Then the probability distribu- 
tion of the ground states in the five subspaces, as well as 
its nature, will be different for taking different parameter 
values. 
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Figure 3: (Color online) Occupation probabilities of the 
ground state in the five subspaces. The inset shows the prob- 
ability distribution of excitations among states with purely 
photonic, purely atomic, and mixed of them. The hopping 
strength between the two cavities is weak for A — O.lg, the 
detuning between atom and field is A = 0, and the atom- atom 
coupling strength is J = g. 



III. THE NATURE OF THE GROUND STATE 



Figure 1: Energy levels for the two-cavity system in the 
absence of hopping (A=0), AEi (i — 1,2,3,4) is the energy 
gap between the adjacent energy levels. 
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In the weak cavity-cavity coupling limit, we can un- 
derstand the nature of the ground state by considering 
the effects of parameters in the effective Hamiltonian 
in Eq.(4). In fact, the energy gap AEi (i = 1,2,3,4) 
between the five subspaces not only dependents on the 
atom- field coupling strength g and detuning A, but also 
relies on the atom-atom coupling strength J, which is 
shown in Fig.l. 

where 



AEi = V(A + J) 2 



i[V(A + J) 2 + 16^ + (A + J)], 



AE 2 

AE 3 
AE 4 



-[(A + J) + V(A + J) 2 + 160 2 ], 
i[V(A + J) 2 + 16^ 2 -(A + J)], 



V(A + J) 2 + Sg 2 

i[V(A + J) 2 + 16s 2 -(A + J)]. (3.1) 



Figure 2: (Color online) Occupation probabilities of the 
ground state in the five subspaces. The inset shows the prob- 
ability distribution of excitations among states with purely 
photonic, purely atomic, and mixed of them. The hopping 
strength between the two cavities is weak for A = 0.1 g, the 
detuning between the atom and the field is A = 0, and the 
atom- atom coupling strength is J = OAg. 



To illustrate the behavior of the ground state operating 
in the dipole-dipole interaction, we first pay our attention 
to the resonant condition A = 0, and three cases for 
J <C g, J ~ g, and J ^> g will be discussed in the 
following. 

At J <C g, there is a large energy gap between 
and |02) due to the photon blockade effect, for which the 
presence of one photon in the cavity blocks the enter- 
ing of the subsequent photons. On this condition, the 
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Figure 4: (Color online) (Occupation probabilities of the 
ground state in the five subspaces. The inset shows the prob- 
ability distribution of excitations among states with purely 
photonic, purely atomic, and mixed of them. The hopping 
strength between the two cavities is weak for A = O.lg, the 
detuning between atom and field is A = 0, and the atom- atom 
coupling strength is J = lOg. 



Figure 6: (Color online) Occupation probabilities of the 
ground state in the five subspaces. The inset shows the prob- 
ability distribution of excitations among states with purely 
photonic, purely atomic, and mixed of them. For large pos- 
itive atom-field detuning, the effects of the dipole coupling 
strength can be ignored. Here, we choose A = O.lg, A = 10 g 
and J = g. 
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Figure 5: The energy gap AEi versus dipole coupling 
strength J at A = 0. The hopping strength between the 
two cavities is weak for A = 0.1 g. 



ground state of the system is approximately (g) \1^ ), 
as shown in Fig. 2. Moreover, is nearly the maximal 
entanglement state of atom and field. For this state, only 
one excitation in each cavity with almost equal probabil- 
ities for atomic and field excitations. This can also be 
confirmed from the inset of Fig. 2. Thus the ground state 
is a polaritonic insulatorlike state, which is analogous to 
the Mott insulator state in the Bose-Hubbard model. 

Compared to Fig. 2, Figure 3 indicates a more different 
behavior for J & g. Besides |0i), the subspace |02) is 
also occupied for the ground state of the system. From 
its inset, we find that both the atom and the field are 




Figure 7: (Color online) Occupation probabilities of the 
ground state in the five subspaces. The inset shows the prob- 
ability distribution of excitations among states with purely 
photonic, purely atomic, and mixed of them. For large nega- 
tive atom-field detuning, two cases of weak and strong dipole 
coupling strength are shown. Here, we choose A = O.lg and 
A = —lOg. Other parameter value is J — O.lg. 



excited, indicating a polaritonic superfluidlike state. 
When J ^> g, the ground state occupies the subspaces 
and |02) with almost same probabilities, as repre- 
sented in Fig. 4. However, only the photons are excited 
in this state. This is because, in this limit, \n~) « —\gn), 
so the ground state is a delocalized photon state in na- 
ture. The state of this form is a photonic superfluidlike 
state. 

The results can be understood easily. Similar to the 
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Figure 8: (Color online) Occupation probabilities of the 
ground state in the five subspaces. The inset shows the prob- 
ability distribution of excitations among states with purely 
photonic, purely atomic, and mixed of them. For large nega- 
tive atom-field detuning, two cases of weak and strong dipole 
coupling strength are shown. Here, we choose A = O.lg and 
A = —lOg. Other parameter value is J = 10 g. 
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Figure 9: (Color online) For weak hopping strength A = O.lg, 
ANi is plotted versus atom-field detuning A and dipole cou- 
pling strength J in the ground state of the coupled two-site 
and two- excitation system. 



effect of detuning in Ref.[29[, at A = 0, the energy gap of 
AEi is a monotonic decreasing function of J, see Fig. 5. 
With the increase of J, the photon blockade effect is de- 
stroyed, leading to an increase of the occupied probability 
of |02). When J & 10g, AEi is almost zero, so the sub- 
spaces and |02) nearly degenerate. Oppositely, AE2 
is a monotonic increasing function of J with a large initial 
value. Thus, the subspace \4>s) can not be occupied. 

Next, we discuss the case of large positive detuning. 
At A = 0, the nature of the ground state changes with 
different values of J. However, when A ^> g, no matter 
what value of J, the energy gap of AEi is always zero. 
While AE2 is a monotonic increasing function of A + J. 
Even J = 0, the energy gap AE2 is very large. There- 
fore, the ground state has equal occupation probabilities 
in the subspaces |0i) and Ifa), but zero probability in 
other subspaces. More interestingly, in the large posi- 
tive detuning limit, the superposition coefficients of the 
dressed states in Eq.(5) have particular values, sin ^-^0 
and cos^r^ 1. Then, \n~) « —\g)\n). Consequently, not 
only the constitution of the ground state but also its na- 
ture is fixed. In Fig. 6, the inset shows that the excitation 
is photonic rather than atomic. This is because the en- 
ergy of the atoms is larger than that of the photons when 
A > 0. Analogous to the situation shown in Fig. 4, the 
ground state indicates photonic superfluidlike nature. 

At A < 0, the energy of atoms is smaller than that of 
photons. In the limit of —A ^> g, J, we find AE\ ~ |A|, 
the ground state is approximate (g) 1 1^" ) . In this con- 
dition, sin%^ 1 and cos ^ 0. Thus, |1") « |e)|0), 
The excitation is almost atomic rather than photonic, 
standing for an atomic insulatorlike state. The result is 
illustrated in Fig. 7. However, when J approaches to the 



value of —A, the value of | A + J\ decreases. That is, the 
energy difference between atom and photon is less and 
less. As a result, the photonic and atomic excitations co- 
exist, as shown in Fig. 8, corresponding to the polaritonic 
insulatorlike state. 



IV. PHASE DIAGRAMS 

The phase diagrams of these states can be distin- 
guished using the corresponding "order parameters" . In 
superfluid states, the excitations in each cavity is uncer- 
tain, resulting in a nonzero variance of the total exci- 
tation number AN\. Oppositely, in the insulator state 
the number of excitations per cavity is constant and thus 
has zero variance. However, the insulator state may be 
either atomic or polaritonic and the superfluid state may 
be photonic or polaritonic in nature. To determine the al- 
lowed types of particles involved in the state, the atomic 
excitation number variance ANai should be taken as the 
"order parameter". AA^u = is corresponding to the 
atomic insulator state or the photonic superfluid state, 
while ANai > revealing the polaritonic nature. Only 
when AA/iAA^u > 0, it shows the polaritonic superfluid 
characteristic for the ground state of the system. 

To begin with, we should distinguish the insulator and 
superfluid areas in the phase diagram, which are deter- 
mined by the total excitation number AN±. With Eq.(4), 
we can give out the expression of ANi directly. 



AiVi 



a\a\ 



{nD - (my 



(4.1) 



In Fig. 9, ANi is plotted under a wide range values of the 



6 




the second "order parameters" ANai, where 



Figure 10: (Color online) For weak hopping strength A — 
O.lg, ANai is plotted versus atom- field detuning A and dipole 
coupling strength J in the ground state of the coupled two-site 
and two- excitation system. 
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Figure 11: (Color online) For weak hopping strength A — 
O.lg, ANi • ANai is plotted versus atom-field detuning A and 
dipole coupling strength J in the ground state of the coupled 
two-site and two- excitation system. 



atom-cavity detunings A and the dipole-dipole interac- 
tion strength J. It is apparent that the phase diagram is 
divided into two sections. The insulator region is under 
the boundary where A < while above the boundary 
is the superfluid region. There is also an area which is 
symmetric to the insulator area, where the ATVi with a 
maximum value 0.5, indicating a most evident superflu- 
idity. 

To find out the polaritonic area, we should introduce 



Nai 
AN A1 



a\a\ 
(N 2 A1 ) 



(j\(Tl, 



(N A1 ) 2 . 



(4.2) 



In Fig. 10, we find that the polaritonic area approximately 
spreads at both sides of the line J = —A. The more 
closer to the line the more obvious the polaritonic na- 
ture. In fact, when J = — A, AEi = (2\/2 — 2)g, it is the 
conditions for photon blockade effect obviously. Then 
the ground state only occupies the subspace \<fri), with 



sin ^j- = cos ^ = ^, standing for a maximal entangle- 



ment state of the atom and the field. So the ground 
state indicates polaritonic insulatorlike nature. 

So far we may guess that on the bases of what are 
embodied in Fig. 9 and Fig. 10, there must be some over- 
lapped areas in the two figures. In Fig.ll, the product 
of ANi and ANai is shown as a contour plot. We can 
clearly identify that there is an area not only in the su- 
perfluid region in Fig. 9, but also in the polaritonic area in 
Fig. 10, representing polaritonic superfluidlike characters. 

Then it is obvious that the phase space in Fig. 9 is 
divided into four sections, from left bottom to top right 
corner, in order of atomic insulatorlike region, polaritonic 
insulatorlike region, polaritonic superfluidlike region and 
photonic superfluidlike region. 



V. CONCLUSIONS 

In summary, we have investigated the QPT of a system 
composed of two coupled cavities, each containing a pair 
of two-level atoms with dipole-dipole interaction. In the 
conditions of fixed cavity-cavity interaction and atom- 
cavity coupling strength, the nature of the ground state 
is dependent on the constituents of the dressed states 
in each cavity and the occupation probabilities of the 
ground state in the five subspaces. Moreover, both of 
them attribute to the dipole-dipole interaction strength 
between the localized atoms and the atom-field detuning 
in each cavity. By choosing three different order param- 
eters, we found that the ground state of the system rep- 
resented more richer behaviors than the Bose-Hubbard 
model. Four types of states are revealed, which divide 
the phase space into four regions. They are the atomic 
insulatorlike state, the polaritonic insulatorlike state, the 
polaritonic superfluidlike state and the photonic super- 
fluidlike state. In the scope of parameter values we taken 
in this paper, the insulator or superfluid phases is deter- 
mined by the combinative effect of A and J, that is the 
value of A + J. Small negative values of it is in favour 
of polaritonic insulatorlike states while for small positive 
value of it embodies polaritonic superfluidlike state. The 
more larger negative value of A + J, the more obvious 
of atomic insulator like nature, and oppositely it shows 
photonic superfluidlike nature. 
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